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Abstract 


The electrical properties of extracellular space around neurons are important to un¬ 
derstand the genesis of extracellular potentials, as well as for localizing neuronal ac¬ 
tivity from extracellular recordings. However, the exact nature of these extracellu¬ 
lar properties is still uncertain. We introduce a method to measure the impedance 
of the tissue, and which preserves the intact cell-medium interface, using whole-cell 
patch-clamp recordings in vivo and in vitro. We find that neural tissue has marked 
non-ohmic and frequency-filtering properties, which are not consistent with a resistive 
(ohmic) medium, as often assumed. The amplitude and phase profiles of the measured 
impedance are consistent with the contribution of ionic diffusion. We also show that 
the impact of such frequency-filtering properties is possibly important on the genesis of 
local field potentials, as well as on the cable properties of neurons. The present results 
show non-ohmic properties of the extracellular medium around neurons, and suggest 
that source estimation methods, as well as the cable properties of neurons, which all 
assume ohmic extracellular medium, may need to be re-evaluated. 


Introduction 

The genesis and propagation of electric signals in brain tissue depend on its electric 
properties, which can be simply resistive (ohmic) or more complex, such as capaci¬ 
tive, polarizable or diffusive. The exact nature of these electric properties is important, 
because non-resistive media will necessarily impose frequency-filtering properties to 
electric signals [1, 2], and therefore will influence any source localization. These elec¬ 
tric properties were measured using metal electrodes, which provided measurements 
suggesting that the brain tissue is essentially resistive [3, 4, 5]. However, the electrical 
behavior of tissue and electrodes can be easily confused [6]; efforts in the direction of 
distinguishing or separating them abound ([7, 8, 9]). Another experimental approach 
using very low-impedance probes revealed marked frequency dependence of conduc¬ 
tivity and permittivity [10, 11]. Indirect evidence for non-resistive media were also 
obtained [12,13,14], and also indicated a marked frequency dependence. 

To explain these discrepancies, we hypothesize that the apparently contradictory 
results are due to the fact that different measurement methods were used. The use 
of metal electrodes represents a non-physiological interface for interacting with the 
surrounding tissue, while in reality, neurons interact with the extracellular medium 
by exchanging ions through membrane ion channels and pumps. To respect as much 
as possible these natural conditions, we have measured the tissue impedance using a 
neuron as a current source, thereby respecting the natural interface. This intracellular 
measurement provides a measurement of a global cell impedance, which contains the 
membrane impedance and the impedance of the extracellular medium. This global 
intracellularly-measured impedance can also be defined as the impedance as seen by 
the cell. 

Figure 1 illustrates this concept and the recording setup needed to record this global 
impedance intracellularly, in vitro (Fig. 1A) or in vivo (Fig. IB). Fig. 1C gives two circuit 
configurations for this system, emphasizing three impedances: Z r is the impedance of 
the intracellular medium (cytoplasm), Z e is the extracellular impedance, and Z RC is 
the membrane impedance, represented by a simple resistance-capacitance (RC) circuit 
(left), or a more complex circuit including dendritic compartments all described by dif- 
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ferent RC circuits (right). An intracellular electrode will measure a global combination 
of these impedances. In the following of the text, we will call this global impedance the 
"global intracellular impedance". 

A central point of our study is that this global intracellular impedance is different 
than the electrical impedances measured by metal electrodes, which we refer here as 
"metal-electrode impedance". We will investigate whether the intracellularly-measured 
impedance reveals more complex electrical properties than with metal electrodes, which 
could possibly explain the discrepancies described above. The global intracellular 
impedance provides not only a realistic estimate of the electrical properties the ex¬ 
tracellular medium, but it is also closer to the natural conditions, and could be a useful 
physical parameter to determine a more precise source localization of cerebral electric 
signals, and to model the propagation of electrical signals in the extracellular space or 
in dendritic trees, as we illustrate here. 


Materials and Methods 

Animals 

Maintenance, surgery and all experiments were performed in accordance with the local 
animal welfare committees (Center for Interdisciplinary Research in Biology and EU 
guidelines, directive 2010/63/EU, and regional ethics committee "Ile-de-France Sud" 
(Certificate 05-003)). Every precaution was taken to minimize stress and the number 
of animals used in each series of experiments. Animals were housed in standard 12 
hours light/dark cycles and food and water were available ad libitum. 


In vitro electrophysiology 

Brain slice preparation. 

Visual cortex. 300 pm-thick coronal brain slices of juvenile mice (P 12-16 Swiss mice 
bred in the CNRS Animal Care facility, Gif-sur-Yvette ; French Agriculture Ministry 
Authorization: B91-272-105) were obtained with a Leica VT 1200S microtome (Leica 
Biosystems, Wetzlar, Germany). Slices were prepared at 4 °C in the following medium 
(in mM): choline chloride 110, KC12.55, NaH 2 P0 4 1.65, NaHCCh 25, dextrose 20, CaCl 2 
0.5, MgC 1 9 7. Slicing started 2 mm from the posterior limit of olfactory bulb, and ended 
3.9mm further. Before recordings, slices were incubated at 34°Cin artificial cerebro¬ 
spinal fluid (ACSF) of the following composition (in mM): NaCl 126, KC1 3, NaHC0 3 
26, NaH 2 P0 4 1.25, myo-inositol 3, sodium pyruvate 2, L-ascorbate de sodium 0.4 ; 
dextrose, 10. The slicing and recording solution was bubbled with 95% 0 2 and 5% 
C0 2 , for a final pH of 7.4. 

Dorsal striatum. Horizontal brain slices with a thickness of 330 pm were prepared from 
rats (P23-30 OFA rats (Charles River, L'Arbresle, France), using a vibrating blade mi¬ 
crotome (VT1200S, Feica Microsystems, Nussloch, Germany). Brains were sliced in a 
95% 0 2 / 5% C0 2 -bubbled, ice-cold cutting solution containing (in mM): NaCl 125, 
KC1 2.5, glucose 25, NaHC0 3 25, NaH 2 P0 4 1.25, CaCl 2 , MgCl 2 1, pyruvic acid 1, and 
then transferred into the same solution at 34 °C. 

Electrophysiological recordings. 

Patch-clamp recordings in pyramidal cells of visual cortex from mice were performed 
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as followed. Slices were superfused at 2 mL/min with the same ACSF solution that 
was used for incubation. Bath temperature of the submerged chamber was main¬ 
tained at 34 °C using a TC-344B temperature controller (Warner Instruments Company, 
Hamden, CT, USA). Neurons in slices of the mouse visual cortex (Pi 2 _i 6 ) were identi¬ 
fied with an upright microscope (Axioscope FS, Zeiss, Germany), an infrared camera 
(C750011, Hamamatsu Photonics, Japan) and an infrared filter. Patch-clamp in the 
whole-cell current clamp configuration in layer V pyramidal neurons was performed 
simultaneously with an extracellular recording using a 3 MU patch pipette (Fig. 1A). 
The latter was located within a close vicinity (~ 30 //m) of the patched cell. All results 
were indifferent to the distance between the reference electrode and the patched cells, 
as potential variations on the extracellular electrode did not excede 1% of the vari¬ 
ations of the intracellular potential. Borosilicate pipettes (1B150F4, World Precision 
Instruments, Inc., Sarasota, Florida, U.S.) of 5-8 MU impedance were used for whole¬ 
cell recordings and contained the intracellular solution (in mM): Hepes 10, EGTA 1, 
K-gluconate 135, MgCl 2 5, CaCl 2 0.1, with osmolarity of 308 mOsm and a pH of 7.3. 
Serial resistance was not compensated for. Recordings were performed with a Multi¬ 
clamp 700B amplifier (Axon Instruments Inc., California, U.S.), filtered at 10 kHz with 
a built-in Bessel filter, and sampled at 25 kHz. Data acquisition and stimulation were 
performed with a National Instruments BNC 2090 A card, and the software ELPFIY (G. 
Sadoc, CNRS, UNIC, France). 

Patch-clamp recordings in medium-sized spiny neurons of dorsal striatum from 
rats were performed as previously described [15]. Briefly, borosilicate glass pipettes of 
6-8M impedance contained for whole-cell recordings (in mM): K-gluconate 105, KC1 
30, HEPES 10, phosphocreatine 10, ATP-Mg 4, GTP-Na 0.3, EGTA 0.3 (adjusted to 
pH 7.35 with KOH). The composition of the extracellular solution was (mM): NaCl 
125, KC1 2.5, glucose 25, NaHC0 3 25, NaH 3 P04 1.25, CaCl 2 2, MgCl 2 1, pyruvic acid 
10, bubbled with 95% 0 2 and 5% C0 2 . Signals were amplified using EPC10-3 am¬ 
plifiers (HEKA Elektronik, Lambrecht, Germany). All recordings were performed at 
34 °C using a temperature control system (Bath-controller V, Luigs & Neumann, Ratin- 
gen, Germany) and slices were continuously superfused at 2 ml/min with the extracel¬ 
lular solution. Slices were visualized on an Olympus BX51WI microscope (Olympus, 
Rungis, France) using a 4x/0.13 objective for the placement of the stimulating electrode 
and a 40x/0.80 water-immersion objective for localizing cells for whole-cell recordings. 
During the experiment, individual striatal and cortical neurons were identified using 
infrared-differential interference contrast microscopy with a CCD camera (Optronis 
VX45; Kehl, Germany). Serial resistance was not compensated for. Current-clamp 
recordings were filtered at 2.9 kHz and sampled at 16.7 kHz using the Patchmaster 
v2x73 program (HEKA Elektronik). Stimuli in current-clamp mode underwent a high 
cut 10kHz filter before being applied. Recordings were made with EPC 10-3 amplifiers 
(HEKA Elektronik; Lambrecht, Germany) with a very high input impedance (1 T U) 
to ensure there was no appreciable signal distortion imposed by the high impedance 
electrode (Nelson et al., 2008). Sinusoidal stimuli were then introduced in whole-cell 
patch-clamp to the patched cell. 

This configuration enables estimating the extracellular impedance, according to the 
circuit displayed in Fig. 1C. To this end, a white noise current stimulus was injected 
into the recorded cell and the impedance was calculated based on this current injection 
(see Fig. 2A). 20 to 120 seconds of Gaussian white noise with zero mean and 100 pA 
variance was injected (results were similar for 30 and 200 pA variance). For each cell, 
we injected 15 to 30 times the same sequence of white noise ("frozen noise") and aver¬ 
aged the measured voltages. This enhanced the signal to noise ratio without altering 
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the results. Figure 2A (right panel) shows the very low autocorrelation of the injected 
currents, being thus a good approximation of white noise. To verify that the same re¬ 
sults can be obtained using a different protocol, we also performed slice experiments 
using sinusoid current stimuli, at different frequencies. This set of experiments was 
performed following the methods previously published [14,16]. Namely, sine waves 
of 12 different frequencies were tested, varying approximately evenly on a logarith¬ 
mic scale ranging from 6 Hz to 926 Hz. Up to 500 traces of 100 to 1500 ms in length 
were averaged before recording the data to disk for offline analyses. Longer stimulus 
lengths and more traces were necessary for the low frequency stimuli. The order of 
the presentation of the frequencies was randomized. Stimuli were introduced with the 
patch electrode in current-clamp mode. The injected current amplitudes ranged from 
200 to 300 pA. Before conducting experiments, we verified via control recordings with 
an external signal generator in the bath without a slice that any amplitude changes or 
phase shifts in the recordings across frequencies induced by the amplifier and record¬ 
ing hardware were negligible [16]. 


In vivo electrophysiology 

Surgical preparation of animals. 

Adult rats (P40-90) were placed in a stereotaxic apparatus (Unimecanique, Asnieres, 
France) after anesthesia induction with a 400mg/kg intra-peritoneal injection of chlo¬ 
ral hydrate (Sigma-Aldrich, Saint-Quentin Fallavier, France). A deep anesthesia main¬ 
tenance was ensured by intra-peritoneal infusion on demand of chloral hydrate deliv¬ 
ered with a peristaltic pump set at 60mg/kg/hour turned on one hour after induction. 
Proper depth of anesthesia was assessed regularly by testing the cardiac rhythm, EcoG 
activity, the lack of response of mild hindpaw pinch and the lack of vibrissae move¬ 
ment. The electrocardiogram was monitored throughout the experiment and body 
temperature was maintained at 36.5 °C by a homeothermic blanket. 

Two craniotomies were performed, one for the insertion of a reference electrode in 
the somatosensory cerebral cortex (layer2/3) and one to allow the whole-cell record¬ 
ings in the contralateral cortex (layer 5). For whole-cell recordings, a 2x2 mm cran¬ 
iotomy was made to expose the left posteromedial barrel subfield at the following co¬ 
ordinates: posterior 3.0-3.5 mm from the bregma, lateral 4.0-4.5 mm from the midline. 
The dura was opened and the craniotomy was filled with low melting point paraffine 
after each time lowering a recording pipette. To increase recording stability the cisterna 
magna was drained. 

Electrophysiological recordings 

Borosilicate glass pipettes of 5 — 8 MU impedance for blind whole-cell recordings 
contained (in mM): K-gluconate 105, KC1 30, HEPES 10, phosphocreatine 10, ATP-Mg 
4, GTP-Na 0.3, EGTA 0.3 (adjusted to pH 7.35 with KOH). Signals were amplified using 
an EPC10-plus-2 amplifier (HEKA Elektronik). Series resistance was not compensated 
for. Current-clamp recordings were filtered at 2.5 kHz and sampled at 5 kHz using 
the Patchmaster v2x32 program (HEKA Elektronik). Whole-cell recordings were per¬ 
formed in pyramidal cells of the somatosensory cortex in layer IV /V (depth from the 
dura: 0.8-1.2 mm) (Fig. IB). Recorded cells were identified as pyramidal cells accord¬ 
ing to their characteristic spiking pattern. The reference was a silver wire placed in 
the contralateral hemisphere. Note that for both in vivo and in vitro experiments, the 
reference electrode is passive, just measuring the extracellular voltage, and thus the 
exact nature of this reference is not critical. Accordingly, the same configuration using 
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a silver micro wire gave similar results in vitro (not shown). 

It is important, however, that the reference electrode be placed in the brain tissue, 
so that the estimated extracellular impedance is not influenced by other tissues. Thus, 
as in the in vitro experiments, this intracellular-extracellular configuration enables es¬ 
timating the extracellular impedance. An important difference with in vitro, is that in 
vivo the current can flow more freely in 3 dimensions, and is closer to natural condi¬ 
tions. Another difference is that in vivo, the system is not silent but displays prominent 
spontaneous activity. To limit this contribution, we have used a "frozen noise" proto¬ 
col, where identical sequences of white noise stimuli were injected repetitively, and the 
sequences averaged. 

The main purpose of this experiment was to estimate Z eq , the equivalent impedance 
between the Ag-AgCl electrode and the ground. In the case of a simple, single-compartment 
neuron, it can be formally defined as Z eq — Z, + Z RC + Z e . This was achieved by ap¬ 
plying two protocols of subthreshold current injection: 

1. The "frozen noise" protocol consisting in the same template of 20 seconds of a 
white noise current (repeated 50 times with 2 s intervals between stimulations 
and averaged). Sequences in which spikes were induced were discarded. 

2. Sinusoidal current at fixed frequencies ranged from 6 Hz to 926 Hz (similar to 
those used in in vitro experiments). The order of the presentation of the frequen¬ 
cies was randomized. 

For the frozen noise and sinusoidal stimuli, the injected currents were tuned for 
each neurons to evoke voltage response of magnitude ranged between 2 and 6 mV. 

Note that in some experiments we injected a hyperpolarizing current (of maximum 
amplitude —150 pA) to prevent suprathreshold activity during application of stimuli. 


Analysis 

All analyses were performed using Python (Python Software Foundation, Wolfeboro 
Falls, New Hampshire, USA), the Scipy Stack and Spyder (Pierre Raybaut, The Spyder 
Development Team). 


In vitro and in vivo patch-clamp - sine waves experiments 

For each frequency and current intensity, the recorded voltages and injected intensities 
were averaged and fit with sines using the optimize package included in Scipy. The 
adequation between data and the fitted sine waves was checked by human interven¬ 
tion for every set of data. The voltage and current were represented respectively as 

V(t, f ) = V 0 sin(2nft + 4 > v ) and I(t, f ) = I 0 sin(2nft + </>*)). The impedance for a given 
frequency was thus given by Z(f) = U e *(V -<?;)_ 


In vitro and in vivo patch-clamp - white noise experiments 

Several models were fit to the experiments. First, we used a purely resistive model (Fig. 
1, bottom left) in which the intracellular impedance Z, can be considered zero and the 
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extracellular impedance Z e was a small resistance (R e ). The equivalent impedance is 
thus given by: 


Zeq, i(^) — 




1 + iu)T n 


+ Re 


(1) 


Second, we combined Zi and Z e as a diffusive term. Rather than trying to find an 
elusive general solution for the usual Nemst-Planck equation [17], we used a first- 
order approximation for ionic diffusion. The impedance of an electrolyte showing non- 
negligible ionic diffusion was derived by Warburg [18, 19], and yielded a modulus 
scaling in 1/y/uj and a constant phase. Similar derivations have been performed in 
different symmetries [20, 21, 22], Note that the latter derivations model the impedance 
of ionic accumulations close to the membrane, by using a first-order approximation of 
the electric potential generated by ionic species following Boltzmann distributions. 

This "diffusion" impedance has been observed experimentally (reviewed in ref. [7]), 
and can be modeled in spherical symmetry by two components scaling the modulus 
and phase ( A and B), and a cutoff frequency f w = lo w /2tt. 


Z w (uj) 


A + iB 
1 + Uj/uJw 


( 2 ) 


This leads to the following expression for the equivalent impedance: 


Z e q,2(u) — 


Rr, 


A + iB 


1 + iUJTn 


1 + Ji— 

1 ' LOW 


(3) 


To take into account the fact that some of the current can flow through the dendrites 
of the cell, we define a dendritic input impedance Z d (see Fig. ID, right), namely the 
impedance of the dendritic tree seen by currents leaving the soma. These currents 
will flow downwards the gradient potential from the intracellular potential ( Vi ntra ) to 
the reference ( V re f = 0). Thus, Z d is defined by Vintra/ ^ where i g d is the generalized 
axial current in the dendrite at the level of the soma, and V intra is the intracellular 
potential at the soma. Note that we consider here Vi ntra as not necessarily equal to 
the transmembrane potential (V m ) because we take V e , the potential of the extracellular 
medium, into account. We then consider separately the resistance and capacitance of 
the soma ( R S ,C S ) and the impedance of the dendrite. The equivalent impedance is 
then: Z eq)3 _ A (oj) = z^+z + Z t + Z e , where Z s (u) = 1+ ^; r is the impedance of the soma, 
Z eq>3 and Z eqA are the impedances for resistive and diffusive media respectively. A 
description of these different models can be found in the next section; parameters for 
each of these models are listed in Table 1. 


Fitting models with and without dendrites 

Two types of models were fit to the experimental measurements, as illustrated in Fig. 1C: 
a single-compartment model, and a model including a dendritic segment. Dendritic 
filtering has been proposed to explain the frequency-dependence of LFPs [23, 24]; 
thus, current flowing in the dendrites could be involved in shaping the measured 
impedance, which would in turn influence the frequency-dependence of local field 
potentials. We tested this possible influence by considering models that include an 
equivalent dendritic compartment, which has been shown to be electrically equivalent 
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to a full dendritic tree [25,26], leading to a "ball-and-stick" type model (see right circuit 
in Fig. 1C). 

In order to fit the models to experimental data, several traditional fitting methods 
were tested (e.g. Newton-Gauss, Levenberg-Marquardt, conjugate gradient, simplex...) 
but were plagued with three main problems: a long computation time (with 4 param¬ 
eters or more), a tendency to get trapped in local minima, and an extreme sensitivity 
to the first estimation of model parameters. We thus developed a probabilistic, non¬ 
iterative method that had none of these problems. 

We proceeded as follows: 1. For each parameter, we defined a domain of accept¬ 
able values, keeping it very large (too much restriction on a parameter is a symptom 
of analytical bias). 2. We drew random sets of parameters (between 500 and 5000), 
e.g Pj = (R m , C m , A, B , fw-), and computed the theoretical impedance spectra they 
predicted with a given model. 3. We computed the squared error Ej(Pj) between the 
theoretical and measured impedances. The error was computed as the sum of squared 
errors on real and imaginary parts.. 4. The couple (Pj, Ej(Pj)) that had the smallest 
errors over all tries was kept as the best fit. 

Random drawing removed the sensitivity to local minima and initial parameters 
that is intrinsic to traditional iterative methods. It allowed a thorough exploration of 
the parameter space, and with high reliability and acceptable efficiency. Empirically, 
this method was found to be much faster than a systematic exploration, but as reliable. 


Modeling the contribution of dendrites 


To model the impedance of the cell including the contribution of dendrites, we use the 
generalized cable formalism (FO model in ref. [28]), which reads: 


2 d 2 V m (x,uj) 
dx 2 


K 2 V m (x,U ) 


(4) 


where 


for a cylindric compartment. Here, the quantity z, is an equivalent impedance, which 
depends on the model considered, z, = r, + r e for the standard (resistive) cable model, 
defined from the intracellular and extracellular resistivities, r, and r e/ respectively. In 
the case of a frequency dependent impedance, z, is more complex and is given by 


1 + iLUTn 


(5) 


Jm) 

Zi = Zi/[1 + ^—(1 + iuJT m )] , 


where z, and z e are the intracellular (cytoplasm) and extracellular impedance densities, 
respectively, r m is the membrane resistivity and r m is the membrane time constant. 
The estimation of these parameters from the experimental measurements is given in 
Supplementary material (Appendix 1). 

When including a dendritic segment, the equivalent impedance (circuit shown in 
Fig. 1C, right) is given by: 



2 (E s + Z e )Zd 

Z s + Z e + Zd 


( 6 ) 


where Z s = 1+ ^ t is the impedance of the somatic membrane, Z%( u) and z e q (u) 
correspond to resistive and diffusive media, respectively. Note that in these models. 
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we have considered Zj & 0 (see Fig. 1C) because the cytoplasm impedance of the soma 
is negligible compared to the membrane impedance. 

If i d is the current flowing in the dendritic tree, the dendritic impedance (as seen by 
the soma) is: 


Taking into account V intra = 


Z d = 


v ; 


intra 


Vintra V m 


9 V 

d v ™ 

V m + V e , we obtain: 


( 7 ) 


14 - 


intra 


Vrr. 


1 + ^L 
V m 


1 + 


Ze 

Z s 


because the conservation law for the generalized current implies V e = Z e (i 9 — i 9 d ) 
and V m = Z s (i 9 — i 9 d ). Note that these equalities would not make sense with the free- 
charge current, because the variations of V m may imply charge accumulation around 
the membrane (dendrite and soma), and thus there is no guarantee of conservation of 
the free-charge current. 


The second part of the fraction represents the input impedance of the dendrite Z, 
which is given by: 

Frn. ^ 


-g- = Z in = — coth(K\ld ) 




where K\ is the cable parameter of the dendrite. Thus: 


Z g \ Ajq 


Zd= 1 + -Z-) —coth(K\.l d ) 




( 8 ) 


( 9 ) 


Note that the values of parameters (k\ and zi m ’ ) in the models considered above 
correspond to an open-circuit configuration which corresponds to the present experi¬ 
ments; the cable equation for the open-circuit configuration, with an arbitrarily com¬ 
plex extracellular medium, was given previously [28]. 

The intermediate formulas and variables for each model are listed in Table 2. 


Statistics on population data 


Different models call for different sets of parameters. For example, the membrane re¬ 
sistance and capacitance of a resistive model are similar but not identical to their coun¬ 
terparts in a model that features a Warburg impedance (e.g., some of the frequency¬ 
filtering properties can come from this supplementary impedance). Thus, for a single 
neuron, we allowed membrane resistance and capacitance to differ across models. 

In order to determine which model was best to describe experimental data, we used 
the following classical procedure. 

1. For each cell and each model, we computed the residual sum of squares (RSS) 
between the experimental curves {y exp ) and theoretical curves (;?/ //,)'■ 

RSS ce u ^ ^ iijexp Ijl.h ) 

frequencies 

For each cell, we normalized this distance by the distance obtained by the best fit. The 
distance between experimental and theoretical curves for a given model was thus: 


RSS rnodr i 


V RSS 


cells 
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2. We want to compare the RSS across models. A raw comparison would be unfair, 
as a model with more parameters is more capable of fitting a given data set. Chosing 
for reference the diffusive model with dendrites (DD), we thus formed for every other 
model M the null hypothesis: "The model M explains the observed curves. If the model 
DD has smaller RSS, it is only because it has more parameters than the model M." 

3. We chose a < 0.05 as threshold for rejecting H 0 . 

4. We used the extra sum of squares F-test, which is able to account for the discrep¬ 
ancy in degrees of freedom across models. We computed the parameter F, following 
an F-distribution under H 0 : 

RSSm - RSS dd DF dd 
RSS dd DF m - DF dd 


where DF is the degrees of freedom (number of cells - number of parameters) of a 
given model. 


5. The F cumulative distribution function (fcdf) allows us to compute the likelihood 
of H q (Fig. 5): 

P = 1 - fcdf(F) 


Ifp < a, we reject the null hypothesis: the diffusive model with dendrites is sig¬ 
nificantly better than the other model, and this can not be explained by the surplus of 
parameters alone. 


Computational models using the measured impedances 


Two types of models were used to test possible consequences of the measurements. 
First, we used a model of the genesis of the extracellular potential. To this end, a current 
waveform corresponding to the total membrane current generated lAp(t ) by an action 
potential was computed from the Hodgkin-Huxley model in the NEURON simulation 
environment [27], This current waveform was used as a current source to calculate the 
extracellular potential, using a formalism that is valid for any extracellular impedance. 
We calculated the extracellular potential by using the impedance measurements made 
in the present paper. According to 2, we have: 


Z e (uj) 


A + iB 

1 + i U / 0 J\y 


( 10 ) 


where A = 151 x 10 6 O, B = 2.54 x 10 6 U and % = 335 rad.s 1 for a distance of a few 
fim. The extracellular potential V ( t ) was calculated using the convolution: 


V(t) 



r)/ A p(r)dr , 


( 11 ) 


where Z e (t) is the inverse Fourier transform of Z e ( uj). 

Second, we simulated a ball-and-stick model using the generalized cable equa¬ 
tion [28] (see Eq. 5). 

A zero-mean Gaussian white-noise current waveform was injected into the den¬ 
drite, and the generalized cable was used to compute the membrane potential in den¬ 
drites and in the soma (see details of the methods in ref. [28]). 
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Results 


We start by outlining the measurement paradigm and the notion of global intracellu¬ 
lar impedance, then we successively describe the results obtained in vitro and in vivo. 
Finally, we illustrate consequences of these findings on two fundamental properties, 
the genesis of extracellular potentials, and the voltage attenuation along neuronal den¬ 
drites. 


The global intracellular impedance 


Here, we explore the hypothesis that the extracellular impedance is fundamentally 
different if measured in natural conditions where the interface between the neuronal 
membrane and extracellular medium is respected, compared to metal electrodes, which 
rely on an artificial metal-medium interface. In natural conditions, the neuronal mem¬ 
brane's interface with the medium involves the opening/closing of ion channels, ionic 
concentration differences and ionic diffusion, whereas metal electrodes involve a dif¬ 
ferent type of ion exchange with the medium, which consists of a chemical reaction 
between the metal and the ions in the medium. To measure the impedance in natural 
conditions, it is therefore necessary to use an intact neuron as the interface with the 
medium, to respect the correct ionic exchange conditions. To do this, we performed 
whole-cell patch-clamp recordings, using neurons as natural current sources in the 
surrounding medium. 

This measurement paradigm is illustrated in Fig. 1A-B and consists of a whole-cell 
patch-clamp recording coupled to a micropipette measuring the potential in the extra¬ 
cellular medium close to the recorded neuron, in vitro (Fig. 1 A) or in vivo (Fig. IB). In all 
cases, the recorded neuron is driven by current injection and serves as a natural current 
source in the medium. In this configuration, relating the signals of the two electrodes 
gives a direct access to the extracellular impedance, as shown by the equivalent circuits 
(Fig. 1C). 

According to this equivalent circuit, we have: 


Vi 


intra 


^eq-i 


;9 


( 12 ) 


where i 9 is the generalized current injected by the patch-clamp electrode. The use of 
the generalized current is essential here because it is the only current that is conserved 
if the media have arbitrarily complex impedances [28], and therefore abides to Kirch- 
hoff's current laws in this general case. Note that the current provided by the current 
generator is also a generalized current because capacitive or non-resistive effects are 
negligible in modern generators. 


In contrast, the conservation of the classic "free-charge current" would apply only 
with resistive impedances. The previous equation gives, for the left circuit (single¬ 
compartment cell): 

Z eq {oj) = Zi + -—■—t-F Z e , (13) 

1 + ZOJTm 

and for the right circuit (cell consisting of a soma and an equivalent dendrite): 


ZeqioJ^) 


Zi + 


(Z s + Z^Zd 
Z s + Z e + Zd 


(14) 


where Z eq is the equivalent impedance of each of the two circuits, Z, and Z e are the 
macroscopic impedances of the cytosol and the extracellular medium, respectively. R m 
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is the global input resistance of the cell and r m is the global membrane time constant 
(left circuit); Z s = 1+ ^ r is the impedance of the soma membrane (right circuit), Z d 
is the input impedance of an equivalent dendrite, as seen by the soma, including the 
extracellular medium surrounding it. 

In the standard model, Z, and Z, are usually modeled by a lumped and resistive 
impedance. In the following, Z eq will be called the "global intracellular impedance" 
of the circuit, because in such a recording configuration the neuron acts as a current 
source in the brain tissue. It represents the impedance of the system as seen by the 
intracellular side of the neuron. 

Z e is the impedance of the extracellular medium as seen by the neuron (extracellu¬ 
lar component of the intracellularly-measured global impedance). We will test here if 
the latter impedance is negligible or constant, as usually assumed. We will consider 
resistive and diffusive versions of this impedance and check which one better fits the 
data. 

Z % includes the impedance of the interface between the tip of the electrode and 
the intracellular medium. This interface will be different in whole-cell patch or sharp- 
electrode recording configurations, because of the location and impedance of the pipettes 
themselves. So, the interpretation of the measured impedances may be different in 
sharp-electrode and whole-cell recordings, and we indeed have observed such dif¬ 
ferences (not shown here). In particular, we made sure that the interface of the silver- 
silver chloride electrodes, used for patching and as a reference, does not contribute sig¬ 
nificantly to the equivalent impedance: it is negligible compared to other impedances 
in the circuit and has little frequency dependence between 1 and 1000 Hz (see Fig. 3). 

In order to measure Z eq , we injected a current intracellularly, and measured the 
intracellular potential V in t r a with respect to a reference. Ideally, this reference is a mi¬ 
cropipette in the extracellular medium (V lfp', see scheme in Fig. 1A). In the in vivo ex¬ 
periments, the reference was a silver electrode inserted in the contralateral somatosen¬ 
sory cortex as in Fig. IB. For subthreshold currents, the system can be considered linear 
in frequency: injecting current at an arbitrary frequency yields voltage variations only 
at that frequency (see Appendix 2 in Supplementary Material). This linearity of the 
system is illustrated in Fig. 4. First, for the subthreshold range of V m considered, the 
membrane V-I relation is linear (Fig. 4A). Second, a combination of four sine-wave cur¬ 
rents of fixed frequency generate V m variations that have strictly the same frequencies 
(Fig. 4B), showing that the linear approximation is valid: sine waves appear as sharp 
spectral lines in Fourier space, two orders of magnitude above baseline. Thus, there 
seems to be no significant impact of non-linear membrane ion channels on the record¬ 
ings. Indeed, no action potentials were present in the recordings analyzed here. 

We used two protocols of current injection, either injection of a series of sinusoidal 
currents of different fixed frequencies (6 to 926 Hz), or injection of a broad-band (white 
noise) current, with a flat spectrum between 1 and 10 kHz. In this case, several in¬ 
stances of the same noisy current trace were injected in the cell ("frozen noise" pro¬ 
tocol), which allowed us to take average and increase the signal-to-noise ratio of the 
measurements. This is particularly useful in vivo, to limit the contamination of the 
measurement by spontaneous synaptic activity, which can be very strong in vivo. 
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In vitro measurements of the global intracellular impedance 

Measurements were first performed in vitro by using an experimental setup consist¬ 
ing of a whole-cell patch-clamp recording of a neuron, together with an extracellular 
recording in the nearby tissue in the cortical slice (see Fig. 2A). Using this recording 
configuration, we computed the global intracellular impedance ( Z eq in Eq. 12; see Ma¬ 
terials and Methods) by using either white-noise current injection, or injected sinu¬ 
soidal currents. The results of a representative cell (N = 31) is shown in Fig. 2B-C. 
Both the modulus amplitude and Fourier phase of the impedance are represented. The 
colored curves in Fig. 2B-C show the best fits of different models to the experimental 
data. One can see that the purely resistive model (RC membrane + resistive extracel¬ 
lular medium, blue curves) does not capture the data. We read from equation 1 that 
\Z eq (u)\ scales as 1/c o in the resistive model, which corresponds to a slope of -1, while 
the experimental modulus yields a slope of —0.5±0.1 (Fig. 2D). The resistive model has 
a phase similar to arctan(fcca) with a minimum of about -90 degrees at high frequencies, 
which contrasts with the -50 degrees observed in the data (Fig. 2E). 

So frequency dependence is clearly different from that predicted by the RC-circuit 
membrane model. The best fits of a model taking into account ionic diffusion [22] can 
account significantly better for most of this frequency dependance in different cells 
(green curves; see also Fig. 5). In particular, the 1/a/c o frequency scaling predicted by 
the diffusive model (equation 3), corresponds to the actually observed -0.5 slope of 
the modulus. The phase modulations can also be remarkably well captured by the 
diffusive model. 

We also tested the possible influence of dendrites, by including an equivalent den¬ 
dritic compartment in the circuit (Fig. 1C, right). This addition could not rescue the 
resistive model, which was still unable to match the observations (Fig. 6). We have 
considered variations of dendritic parameters, including very long dendritic segments, 
different axial and leak conductances, and did not see any significant improvement by 
the addition of dendrites. In the diffusive model, taking the dendrites in account only 
enhanced marginally the agreement between experimental and theoretical curves. Sta¬ 
tistical analysis showed that the improvement in quality from the resistive to the dif¬ 
fusive model was significant, and not only due to a higher number of parameters. 
Furthermore, the apparent smaller number of parameters in resistive models can come 
from hidden assumptions, such as homogeneity and low resistivity of the extracellular 
medium. 

These results were replicated in striatal neurons using purely sinusoidal input cur¬ 
rents from 6 to 926 FIz (see Fig. 7, N = 18), thereby confirming that the global shape 
of the global cell impedance does not depend on the stimulation protocol. In addition, 
we tested a capacitive (RC) model of the extracellular space, but this model did not 
account for the modulus and the phase modulations (it was the worst fit of all models 
tested - not shown). 

We also checked whether the quality of seals could affect the global intracellular 
impedance measurements. Indeed, if the cell membrane is bypassed, the impedance 
is not measured anymore through the natural interface of a neuron membrane. The 
average of neurons with good seals (> 1 GU) yields a slope of —0.5 ±0.1 (see fig. 2D); 
in comparison, cells with extremely poor seals (e.g. 200 MU, not included in the data 
shown here) yielded a flatter impedance, with a slope between 0 and -0.3. This can be 
easily explained by replacing Z RC by a resistance in the expression of Z eq . 

Finally, we checked whether part of the observed frequency dependence could be 
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due to the recording pipette. We found that the frequency dependence of the patching 
pipette and silver-silver chloride electrode is negligible in ACSF (Fig. 3). These mea¬ 
surements show that the observed frequency dependence of the impedance cannnot 
be attributed to the silver electrode interface, and probably stems from the properties 
of the extracellular medium. 

Global intracellular impedance in vivo 

In a second set of experiments, the measurements were performed in vivo with whole¬ 
cell patch-clamp recordings of pyramidal cells of the somatosensory cortex, layer V 
(see scheme in Fig. 8A and details in Materials and Methods). Similarly to Fig. 2, the 
modulus and phase of the impedance were estimated by white noise current injection 
(Fig. 8B-C). Although the data display a high degree of noise (due to spontaneous 
synaptic inputs in vivo), they were in qualitative agreement with in vitro results on N = 
18 cells. The resistive model did not capture the modulus amplitude, nor the phase 
of the global intracellular impedance of the neuron. The diffusive model was able 
to capture the essential variations, both in amplitude (modulus) and phase domain. 
Similar to in vitro measurements, the modulus yields a slope of —0.4 ±0.1 (Fig. 8D), 
and a minimum phase around -50 degrees (Fig. 8E), which significantly deviate from 
predictions of a resistive model. 

As in the in vitro experiments, the addition of an equivalent dendritic compartment 
did not improve these differences. The resistive model with dendrites was also unable 
to account for the measurements, while the diffusive model provided acceptable fits to 
the data. 

One must keep in mind that the in vivo measurements were made in the presence 
of low-frequency spontaneous activity, typical of anesthetized states. This "synaptic 
bombardment" probably explains the mismatch of all impedance models at low fre¬ 
quencies in vivo. Such a mismatch was not present in vitro. 


Possible consequences of these measurements 

Finally, to evaluate possible consequences of these measurements, we have considered 
two situations where the extracellular impedance can have strong consequences. A 
first consequence is the fact that the diffusive nature of the medium will necessarily 
impose frequency filtering properties on extracellular potentials, which affects mea¬ 
surements made with extracellular electrodes. To illustrate this point, we simulated 
extracellular potentials generated by a current source corresponding to the total mem¬ 
brane current generated by an action potential (using the FIodgkin-FIuxley model). 
We then calculated the extracellular potential at a distance from this current source, 
using either a resistive model, or a diffusive model (Fig. 9A). Interestingly, one can 
see that the extracellular signature of the spike has a slower time course in diffusive 
conditions. A similar situation was also simulated using subthreshold noisy excitatory 
and inhibitory synaptic activity (Fig. 9B). In both cases, the nature of the medium influ¬ 
ences the shape and propagation of the local field potential (LFP), for both extracellular 
spikes and LFP resulting from synaptic activity. 

A second possible consequence is on the cable properties of neurons. This point 
was illustrated by simulating a ball-and-stick model subject to injection of a noisy cur¬ 
rent waveform in the dendritic cable (Fig. 9C). As shown in Fig. 9D, the attenuation 
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of the voltage along the dendrite can be drastically different in a diffusive medium 
compared to a resistive medium, as noted previously [28]. Including a diffusive extra¬ 
cellular medium reduced voltage attenuation (Fig. 9D, blue curve), but this reduction 
was the strongest when both intracellular and extracellular media were diffusive (red 
cuves). Thus, the nature of the medium will also influence the shape and propagation 
of potentials in dendrites. 


Discussion 

We have provided here the first experimental measurements of the impedance of the 
extracellular medium, in natural conditions, in vitro and in vivo. We found that, not 
only the estimated extracellular impedance is higher compared to traditional metal- 
electrode measurements, but it is also more frequency dependent. The standard model, 
considering the medium as resistive, can account for metal electrode measurements, 
but not for natural impedance measurements. In contrast, we found that a diffusive 
model can account for most measurements, both in modulus amplitude and in phase. 
We also checked whether the inclusion of dendrites could affect these conclusions, but 
it did not qualitatively change these results. 

It is noteworthy that the present measurements are made between a neuron, and 
a reference electrode in the nearby tissue. Therefore, the current presumably flows in 
the entire tissue, and thus, the impedance measured can be considered "macroscopic". 
From the different experiments realized here, we estimate that the impedance is deter¬ 
mined by the region close to the membrane, within distances of the order of hundreds 
of microns in cerebral cortex. 

The apparent inconsistency with the previous metal-electrode measurements can 
be resolved by considering that each kind of electrode has a specific interface and 
impedance, depending on its physical nature [7], Classical impedance measurement 
studies tackle this problem offline with a normalization by a measure in saline [3,10], 
online by removing the effect of the interface using the saturation due to large cur¬ 
rents [5], or by minimizing the interface [29, 30]. In physiological conditions, neurons 
have an electrical interface with the extracellular medium, as a part of their normal 
environment. This interface should therefore not be removed when using neurons to 
evaluate the impedance of the extracellular medium, as it is one of the keys to explain¬ 
ing the electric field produced by an active cell. 

The system presented here deals with the usual problems of electrode recordings 
(see ref. [29]) in unusual ways, which solves some classical issues but raises new in¬ 
terrogations. First, the electrode - or neuron - must be standardized. It is remarkable 
to see that, despite the considerable cell-to-cell variability in size or morphology, we 
obtained here very consistent measurements, with very similar amplitude and phase 
profiles from cell to cell. These measurements can be captured with accuracy with a 
limited number of parameters, most of which are well known ( R mr C m ...). How the 
cellular morphology influences these results, and why this influence seems so small, 
constitute interesting subjects for future extensions of the present work. Second, the 
spatial scale concerned by potential measurement and current injection is a grey area 
in the litterature (see ref. [31]). We believe that the scale of a single neuron may be as 
relevant as the tip of traditional recording electrodes, of arbitrary size and position in 
an inhomogenous medium, which can affect recordings significantly [14]. Third, the 
interface of the electrode and its behavior must be linear and well understood within 
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the measurement range, which we discussed previously. Provided the system is op¬ 
erated with all necessary precautions, linearity is maintained (Fig. 4); the path of the 
injected currents in neuron compartments does not seem to be a crucial matter (Fig. 14) 
and injected currents splitting between different forms (free or bound charges, electric 
flux...) is not a problem within the generalized current formalism. Furthermore, the 
traditional 4-electrode setup is designed to separate voltage recording from possible 
filtering by the interface generated when injecting current [30]. In the system pre¬ 
sented here, the silver-silver chloride wire has a very resistive interface (Fig. 3) and is 
negligible with respect to the main, relevant interface of the recorded neuron. 

A possible explanation for the prominent role of ionic diffusion is that when a neu¬ 
ron acts as a current source, the electric field lines will not, in general, match the com¬ 
plex geometry of the extracellular medium. The trajectory of ions would thus be af¬ 
fected by obstacles such as cells and fibers [14], which would yield local variations in 
ionic concentrations. Ionic diffusion would therefore exert an important force on ions 
in the extracellular medium. A linear approximation of this phenomenon allows one 
to model this contribution by a Warburg-type impedance, scaling as ^=. In addition, 
ionic diffusion is involved in membrane potential changes, and participates as well to 
maintaining the Debye layer surrounding the membrane [32]. Taken together, these 
factors could explain why the present measurements are in such good agreement with 
the diffusive model. 

Despite this agreement, the participation of diffusive phenomena can vary with 
age and experimental conditions. As the brain gets older, the extracellular volume 
fraction shrinks, which could make ionic diffusion even stronger and thus reinforce 
the Warburg component of the natural impedance. Furthermore, in vivo tissue may be 
more confined than in vitro, with a similar result on the importance of ionic diffusion. 
One can thus reasonably expect Z w to be stronger in P30 rats in vivo (N = 18) than in 
P12-16 mice in vitro (N = 31). Indeed, the components A and B of Zw are significantly 
stronger in vivo in P30 rats than in vitro in P12-16 mice (comparing the medians ± 
the standard error of the mean: resp. 96 ± 12 Mil vs 183 ± 16 MQ for A and 6.6 ± 
2.7 MVt vs 28 ± 3.7 MVt for B). Thus, the age of the subject and type of recording need 
to be taken into account when using measured values of the natural impedance. In 
particular, it may lead to an underestimation of Zw in this paper, because we mainly 
focused on in vitro recordings in young mice; our conclusions on the importance of 
ionic diffusion are thus rather conservative. It is noteworthy that in between these two 
sets of observations, the Warburg frequency remains the same : 43 ± 3 Hz in vitro vs 
42 ± 3 Hz in vivo. 

Our results do not disqualify the previous measurements, but are complementary. 
We suggest that for all cases where the current sources are generated by natural condi¬ 
tions (i.e., by neurons), the global intracellular impedance should be used. This is the 
case for example when analyzing the LFP signal, or with Current Source Density (CSD) 
analysis. In cases where a metal electrode is used to inject current, the metal-electrode 
impedance would be relevant, for example, in Deep Brain Stimulation paradigms. 

Note that, although the diffusive model accounts very well for the modulus and 
phase variations of the global intracellular impedance, there exists small deviations, 
in particular at high frequencies. The latter may be due to a number of phenomena, 
including variability in neuron geometry or limitations of the linear approximations 
used here. The existence of "shunts" due to the liquid around the electrodes is also 
not to be excluded. Further studies should be designed to identify the contribution 
of such factors, e.g. pharmacological inactivation of nonlinear channels. Two argu- 
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merits suggest that the present formalism is satisfactory, the strong reproducibility of 
results across 31 recorded neurons in vitro, despite intrinsic biological variability, and 
the coherence between the diffusive model and experimental data. 

The exact boundaries of the domain where these results apply are still to be deter¬ 
mined. For example, in figure 9 we are extrapolating into a nonlinear region to make 
implications about the shape of the action potential. We think this extrapolation is ac¬ 
ceptable because non-linear behavior is mostly happening in the highest frequencies, 
barely overlapping with the LFP frequency range (see also Appendix 2 in Supplemen¬ 
tary Material), but one should be aware of that caveat. 

Finally, using computational models, we illustrated consequences of the medium 
non-resistivity on extracellular and intracellular potentials. A number of fundamen¬ 
tal theoretical equations used in neuroscience, such as CSD analysis [33], or neuronal 
cable equations [25, 26] were originally derived under the assumption that the extra¬ 
cellular medium is resistive. If the medium is non-resistive, these equations are not 
valid anymore and must be generalized. Attempts for such generalizations were pro¬ 
posed recently for CSD analysis [22] and cable equations [28], but they were not con¬ 
strained by measurements. The simulations provided here show that including a dif¬ 
fusive impedance based on the present measurements has significant consequences, 
for both extracellular potentials, and for the electrotonic properties of neurons. The 
shape of the extracellular spike may be affected by the nature of the medium (Fig. 9A), 
assuming that one can extrapolate the present results to the nonlinear region of the 
V m . This shows that the sharpness of the extracellular spike may be influenced by the 
properties of the medium, which constitutes another factor that could complicate the 
identification of neurons from spike shape. The dendritic attenuation is also reduced 
in the presence of a diffusive medium (Fig. 9D-E), as shown previously [28]. Extrapo¬ 
lating these results, it seems that the sources estimated by CSD analysis, could greatly 
be affected by the nature of the extracellular medium, which constitutes a direct exten¬ 
sion of the present study. Similarly, source reconstruction methods from the EEG, are 
also likely to be affected by the nature of the medium, and thus, these methods may 
need to be re-evaluated as well. 
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Tables 


Model type 

Resistive 

Diffusive 

No dendrite 

Rm/ C'm/ Re 

Rm/ C'm/ d, B , f\y 

Dendrite 

Rs, C s/ Re/ I'd 

Rs/ C s , A, B, f W , l d 


Table 1: Parameters used for the four models. l d is the length of the dendritic compart¬ 
ment. R e is the extracellular resistance. 


Model type 

Z i Zi ^ s 

Z e 

(m) 

Ze 


Standard 


negligible 

negligible 

f-(l +iu>r m ) 

I m v 

Diffusive 

n r 

A+iB 

0.5r m 


(1 +i)y/^ ls 

1+^/ iuj/uj w 

27raC m (l+i)v^ 

,( m ) 

r /n[ 1 +-^r^ l+iuiTm \ 


Table 2: Intermediate variables. Zi and Z e are respectively the intracellular and extra¬ 
cellular impedances, Ze is the input resistance of the extracellular medium seen by 
the dendrite (in Vt/m), k\ is the cable parameter of the system. Constants: r, and r e are 
the linear density of resistance in the cytoplasm and extracellular medium (estimated 
to 28 x 10 9 and 18 x 10 9 Vt/m), Zi is the linear density of impedance in the dendrite, 
l s and l,i are the length of the soma and dendrite, a is the diameter of the dendrite, r m 
is the time constant of the membrane (~10 ms), C m is the specific membrane capaci¬ 
tance (1CT 2 F/m 2 ). We also have r m = T rn /2na m C rn , where l/r m is the linear density of 
membrane conductance (S/m). 
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Figure 1: Scheme of the experimental setup for measuring the global intracellular 
impedance. A-B: Placement of electrodes in vitro and in vivo. In each case, a cell is 
recorded in patch-clamp whole-cell configuration using a micropipette, where the ref¬ 
erence electrode is a micropipette located in the extracellular medium at a short dis¬ 
tance from the recording (A), or a silver wire in the contralateral hemisphere (B). C: 
Equivalent circuits corresponding to this experimental setup, in two different configu¬ 
rations: the cell is either considered as a single compartment (left), or with dendrites 
(right), resulting in a slightly more complicated circuit. The membrane is modeled as 
a RC-circuit, where R m and C m are the global membrane resistance and capacitance of 
the cell (left), respectively, or R s and C s are the resistance and capacitance of the so¬ 
matic membrane (right). Z, is the macroscopic intracellular impedance (including the 
electrode-cytosol interface), and Z e is the macroscopic extracellular impedance. In the 
resistive models, Z e is modeled as a simple resistance R e ; in the diffusive models, it is a 
function of three parameters (A, B and k, respectively scaling the modulus, phase and 
cutoff frequency of the diffusive impedance; see Materials and Methods). In the right 
circuit, Z d is the equivalent input impedance of the dendritic tree seen by the soma; 
it is a function of Id, the equivalent length of the dendritic tree. In both circuits, V intra 
represents the intracellular potential, and V)/ p the potential in extracellular space. 
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Figure 2: Global intracellular impedance of cortical neurons recorded in vitro (current 
clamp). A: Left: Scheme of the experimental setup; Middle: example signals; A gaus- 
sian white noise current signal (red, amplitude: ±100 pA) was injected repeatedly in 
patched neurons while recording the intracellular potential (blue). For each 20 seconds 
period, we computed the impedance as the ratio of measured voltage on injected cur¬ 
rent ( Z mes = V/ . ntr " ) and averaged on 30 periods. Acquisition was performed at 20 kHz. 
Right: Autocorrelation of the injected white noise current in Fourier space. B: modu¬ 
lus of Z eq (log scale) represented as a function of frequency (/o;/ I0 (/)), for white noise 
current injection. C: Fourier phase of Z eq , in the same experiment. The different curves 
shows the best fit of two models to the experimental data (blue: resistive, green: diffu¬ 
sive). Parameters of the models: i) Resistive: R m = 200MO, C m = 45 pF, R e = 19 MO; 
ii) Diffusive: R, m = 180 MO, C m = 110 pF, A = 99MO, B = 3.8 MO, f w = 36 Hz. A, B 
and f w are parameters of the diffusive impedance scaling respectively its amplitude, 
phase, and cutoff frequency. D: Distribution of the slopes of Z eq fitted between 20 and 
200 Hz (linear fit, red dashed lines on the figures). E: Coordinates of the minima of the 
Fourier phases for each cell. 


23 


































Modulus (Ohms) 


A 



B 



Figure 3: Impedance of a silver-silver chloride electrode inserted in a typical (~ 3 Mil) 
patch pipette, measured in ACSF. Left: modulus (Ohms), right: phase (degrees). 
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Figure 4: Linearity in temporal and frequency space. A: injection of hyperpolarizing 
and depolarizing pulses in a cortical neuron patched in vivo, around resting membrane 
potential. The V-I relation indicated is roughly linear in this subthreshold range. B: 
Sine waves of current with four different frequencies were injected simultaneously in 
a patched neuron in vivo (top, red curve). The voltage was recorded (top, blue curve). 
The modulus of the Fourier Transform of both signals is shown here, as a function of 
frequency (Hz): current in red (A), voltage in blue (V). The frequencies corresponding 
to the peaks on both signals illustrate the linearity in frequency of the system, here 
tested between 1 and 1000 Hz. 
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Figure 5: Average fitting error of the four models. The figure shows the average good¬ 
ness of fit for 4 different models investigated: Resistive, Resistive with Dendrite, Diffusive, 
and Diffusive with Dendrite. The error bars are standard deviations of fitting errors. The 
testing variable used to compare the quality of these fits takes into account their differ¬ 
ent number of variables. Of all models, the diffusive models give the smallest error in 
all cases tested. The best fit is provided by the diffusive model with dendrite, although 
only marginally better than the diffusive model in a single-compartment model. 
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Figure 6: Dendritic contribution to the global intracellular impedance of a cortical 
neuron recorded in vitro (current clamp). Same arrangement of panels as in fig. 2, 
except that the different curves show the best fit of two models to the experimental 
data (blue: resistive, green: diffusive), both models included an equivalent dendritic 
compartment. Parameters of the models: i) Resistive: R s = 240 MQ, C s = 37 pF, 
R e = 21 MQ, Idend = 390 pm; ii) Diffusive: R s = 150 M Q, C s = 89 pF, A = 130Mf2, 
B = —12 MQ, f w = 30 Hz, l dend = 12 pm. A, B and f w are parameters of the diffusive 
impedance scaling respectively its amplitude, phase, and cutoff frequency. 
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Figure 7: Global intracellular impedance measurements of a striatal neuron recorded 
in vitro (current clamp) and stimulated with sinusoid inputs from 6 to 926 Hz. A: 
sine waves of current (top, red; minimum 25 samples and 20 cycles per sample) were 
injected in patched neurons while recording the intracellular potential (top, in blue; 
16.7 kHz sampling). B: modulus of Z eq (log scale) represented as a function of fre¬ 
quency ( log w (f )), for sine wave current injection (600 pA amplitude). C: Fourier phase 
of Z eq in the same experiment. Parameters of the models: i) Resistive: R m = 120 MO, 
C m = 46 pF, R e = 28 MO; ii) Diffusive: R, m = 130 MO, C m = 59 pF, A = 58MO, 
B = 43 MO, f w = 50 Hz. Representative of N = 18 cells. 
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Figure 8: Global intracellular impedance of cortical neurons recorded in vivo in 
current-clamp. A: Scheme of the recording configuration, with a similar representation 
as in Fig. 2A, but performed in rat somatosensory cortex in vivo (with a contralateral 
reference electrode), using white noise current injection. For these experiments, we 
injected 20 times the same gaussian white noise current trace (top, in red), recorded 
at 20 kHz and averaged the intracellular potential (top, in blue). We calculated the 
impedance seen by the neuron. B: Modulus of the impedance obtained, as a func¬ 
tion of frequency. C: Phase of the impedance. Parameters of the models: i) Resistive: 
R m = 230 MQ, C m = 28 pF, R e = 29 MCI; ii) Diffusive: R m = 100 MVL, C m = 21 pF, 
A = 250 MVL, B = 24 MVL, f w = 20 Hz Representative of N = 18 cells. D: Distribution 
of the slopes of Z eq fitted between 25 and 250 Hz (linear fit, red dashed lines on the 
figures). E: Coordinates of the minima of the Fourier phases for each cell. 
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Figure 9: Consequences of the diffusive nature of the extracellular medium. A. Sim¬ 
ulation of LFP in the extracellular medium (10 //m from the soma) following injection 
of current following a spike waveform. The panel compares the extracellular spike 
obtained for a resistive medium (black), compared to a diffusive medium (red), where 
the filtering is evident. B. Same simulation as in A, but using injection of the combined 
current of noisy excitatory and inhibitory (subthreshold) inputs. In this case, the LFP 
obtained in a resistive and diffusive medium are also differentially filtered. C. Scheme 
of a ball-and-stick neuron model where a noisy current waveform was injected in the 
middle of the dendrite. D. Relative voltage attenuation profile obtained (at 5FIz) when 
the neuron is simulated in resistive (black) or diffusive (red, blue) media. Two diffusive 
configurations were simulated. Diffusive 1: diffusive intracellular and extracellular me¬ 
dia (red curves); Diffusive 2: diffusive extracellular medium with resistive intracellular 
medium (blue curves). Using diffusive media results in a reduced voltage attenuation. 
In all cases, the resistive or diffusive media were simulated using the best fits to the 
impedances measured in vitro. 
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Supplementary material 


A Appendix 1: Integrating the global intracellular impedance 
in models 


In this appendix, we relate the impedance Ze in the generalized cable, to the impedance 
measurements reported in the present paper. In the generalized cable [28], the extra¬ 
cellular impedance was modeled by parameter Ze ■ Starting from the expression of 
the extracellular impedance. 




(m) 

Ze 


(m) 

Ze 


A A.'TT /■? 

soma ^ /[±L soma 


(A.l) 


and condisering a single-compartment model, according to Eq. 13, we can write 


A m ) 


A, 


1 Zp Ag 0ma 


R r 


1 + iOJT n 


(A.2) 


Thus, by assuming a typical somatic membrane area, we can estimate Z e , and thus 
also estimate Ze " ] . The other parameters, R, mr r mr Z eq , can also be estimated from the 
present measurements. 


B Appendix 2: Establishing the linearity of the system 

In this Appendix, we explain how to determine the linearity of the system, in temporal 
and frequency space. 

B.l Linearity in Fourier frequency space 

In the Fourier frequency domain, the experiments show that the ratio AV(u;) / I 9 (u>) 
is a bounded function for variations around the resting membrane potential. In these 
conditions, a sinusoid in current gives a sinusoid voltage with the same frequency, and 
with no additional peak in the spectrum. This is true for relatively small variations (a 
few millivolts), keeping the membrane far away from spike threshold. As shown in 
Fig. ID, a combination of sine-wave currents generates a voltage power spectrum with 
peaks at the same position in frequency, and where no additional peak or harmonics 
appear. We can say that in this case, the membrane potential of the neuron is linear in 
Fourier frequency space. This implies that each component of this system in this space 
is also linear, and in particular, the V-I relation of ion channels in the membrane are lin¬ 
ear, because the membrane capacitance is approximately constant (White, 1970). This 
is an expression of Ohm's law, in which the ion channels are equivalent to a resistor, 
with no voltage-dependent effects (see Section B.2). 

To demonstrate this, we note that the ratio between V and I 9 is a continuous bounded 
function in Fourier frequency space, with the constraint 1/(0) = g(u, 0) = 0 for uj y 0 
(the latter condition means that the neuron is at rest when the transmembrane current 
is zero). In these conditions, we have 

V( U )=g( U ,r>(u)). (B.l) 
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We can develop g in Taylor series relative to the current, because the domain of def¬ 
inition of g is necessarily compact in experimental situations. Consequently, we can 
write: 


= V( W ) - V(0) = A (u , i o)/»( w ) + |0( w ,o)/»» + ... 

= b\(u)I 9 (u) + b 2 (u)I g2 {u) + ■ ■ ■ . 


(B.2) 


The impedance is then given by: 

z M = = &i(w) + h(u)I°(u>) + ■■■ . (B.3) 

We can see that, if the spectrum /(cu) is a discrete Fourier spectrum composed of Dirac 
delta functions, then Z(u) cannot be a bounded function when b n 0 for n > 1. Thus, 
we obtain 

AVu = b 1 (u)I 9 (uj) (B.4) 

when the ratio V(oj)/I 9 (oj) is a bounded function. In other words, the system is nec¬ 
essarily linear in Frequency space because the V-I relation does not depend on the 
current amplitude. Note that this independence is only true in the absence of voltage- 
dependent conductances, so it can apply to the subthreshold range, near the resting 
membrane potential. Such a linear dependence also implies that the position of spec¬ 
tral lines is necessarily identical between V ( u ) and I 9 (uj). 

B.2 Linearity of traditional V-I curves 

We now address the question of whether the V-I relation of ion channels is linear when 
these channels are linear in Fourier frequency space, and vice-versa. 

In general, for a membrane containing ion channels, we have: 

V = /(/) , (B.5) 

with V = f(0)= cst for zero current (resting membrane potential). 

We can approximate V as precise as we want using a polynomial of the current, 
because V is necessarily a continuous function of this variable since the electric field 
is finite (E = — W). This is by virtue of the Stone-Weierstrass theorem [34], which 
states that every continuous function defined over a closed and bounded domain, can 
be approximated as close as we want by a polynomial. Thus, for a given population of 
ion channels, we can write 


A V(t) = V(t) - V(0) = aj + a 2 I 2 + • • • . (B.6) 

If we express I as I(t) = e lult , we obtain 

A V(t) = die iut + a 2 e 2iu}t + ■■■ (B.7) 

such that the Fourier transform of the variations of V around V (0) generally gives a 
spectrum very different from that of the current. Indeed, applying the Fourier trans¬ 
form gives 

AV (c o) = ai 5 (uj' — oj) + a 2 S(u' — 2 u) + • • • (B.8) 
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Thus, it is necessary that Vn > 1 we have a n = 0 if we want that the position of the 
spectral lines of AV (ca) is the same as that of 

Moreover, it is evident that if the function / is linear, then the position of the spectral 
lines of AV(u>) is the same as that of /(cu). 

Thus, the linearity in Fourier frequency space implies linearity of the V-I relation of 
the ion channels activated in the range of V where V = /(/). The linearity in Fourier 
frequency space constitutes a full condition of linearity, because the V-I relation can be 
more complex, for example V = f(u, I). 
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